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Abstract 

The Mean Spherical Approximation for fluids is extended to treat the case 
of dense systems interacting via soft-potentials. The extension takes the 
form of a generalized statement concerning the behavior of the direct correla- 
tion function c(r) and radial distribution function g(r). From a detailed 
analysis that views the hard core portion of a potential as a perturbation 
on the whole, a specific model is proposed which possesses analytic solutions 
for both Cj'i. jmb and Yukawa potentials, in addition to certain other remarkable 
properties. A variational principle for the model leads to a relati/ely 
simple method for obtaining numerical solutions. 

Permanent address; The Nuclear Research Center-Negev, Beer-Sheva, Israel. 
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I. INTRODUCTION 

Th« mean spherical approximaclon (MSA) has been applied^ almo«c exclusively 
to dense liquids wl^ose In ce molecular potentials u(r) possess a hard core of 
range 0. For such potentials the MSA Is specified by the equations 


g(r) - 1 h(r) - 0 , 

r < a , 

(U) 

c(r) + 3u(r) - 0 , 

T > J , 

(lb) 


(9 - 1/kgT) 

'.ogether with the Ornsteln Zemlke (OZ) relation, which for a system of average 
density P Is written 

h(r) - c(r) + p 'dr'hdr-r' l)c(r') . (2) 

Here g(r) Is the radial distribution function and c(r) the direct correlation 

function for which (2) Is the defining relation . 

The mean spherical approximaclon Is not a satisfactory model for low 

density fluids since it normally fails to give the correct value for the second 

vlrlal coefficient. An exception to thi-j however Is the system of hard 

spheres, where the MSA is equivalent to the Percus-Yevick (PY) approximation. 

Another exception Is Che dense Coulomb gas (u(r) r where Che MSA gives 

2 

Che known Debye-Hiickel limit In a first order "Inverse range" expansion: It 

3 

also leads Co results that satisfy the Stllllnger-LoveCC conditions . The 
major Interest In the model, however, can be traced Co the fact that It yields 
Interesting analytic solutions for a fairly wide class of systems providing 
they possess, in their interparticle potentials, the fundamental hard-core 
property as represented in(laX On the othei' hand, the evident disadvantage 
of the model has been an apparent lack of any systematic basis for its extension 
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CO non-hard core pocenclals. The purpose of chla paper La Co provide such a 

basla which is being proposed ntore for Che qualicacive insighc ic gives in Che 

cheory of liquid scruccure chan ica q*iancicacive numerical prediccions. Though 

founded ~n a more general acacemenc oc che MSA, (Co be given in Che nexc 

seccion), its pracCical realizacion is very similar. In parcicular a version 

applicable cc sofc pocencials (che "soft" mean spherical approxictacion) is 

incroduced, and in secCion III is compared Co che modified hypemecced chain 
4 

approach . Some che rmodynamic resulcs are considered 'n seccion IV, and a 
variacional principle for che model is given in seccion V. Finally, 4.'pplica- 
Cions and accendanc procedures are given for Che one-componenC plasma in 
seccion VI. Some conclusions are drawn in seccion VII. 

II. THE GENERALIZED MEAN SPHERICAL APPROXIMATION 

4 

Analysis of a large body of compucer simulacion dace compiled for a 
varieCy of disparace buc mainly sofc incerparcicle pocencials has led co a 
conclusion chac from Che scandpoinc oc iCs Chermodynamic consequences, a 
replacemenC of (lb) by 

'dr g(r)(c(r) + 3u(r)] - 0 (3) 

conscicuces a relacively accurace reformulacion of che approximacion, and 
one which in che dense fluid regime of low compressibilicy is valid for any 
pocencial. The meaning of (3) is mosC easily gauged by recalling Che exacC 
relacion^ 

c(0) - -aOp/dp)^ + 2(3U/N) - p jdr g(r)[c(r) + 3u(r)] (4) 

where p is Che pressure, also sacisfying 
3(3p/ap)-. - 1 - p Tdr c(r) , 


( 5 ) 
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and U/N, the internal energy per particle, is given by 

(3U/N) • ^Epj'g(r) 3u(r)dr . (6) 

(We observe at this point that (4) holds for the one component plasma. For 
this system, hov«ver, the presence of a compensating uniform background 
requires us to replace c(r) in (5) by c(r) + 3u(r) and g(r) in (6) by h(r)0 
The generalized mean spherical approximation (GMSA) is therefore specified 
by the equations 


g(r) - 1 + h(r) - 0, 

(r < 7), 

(la) 

^ dr g(r)(c(r) + 3u(r) • 0, 


(Ic) 


whose solution will determine approximations to c(r) and g(r). For potentials 
with a hard cor<} 7 is Just the hard core diameter. For soft core potentials 
7 can be a paramet^'r (see below). Clearly the model can only be considered 
useful if the exact functions c(r) and g(r) satisfy 

Pj 8(f) [c(r) + 3u(r)] dr « c(0) . (7) 

However, js already noted, condition (7) appears to be well satisfied in 
many dense systems of interest. 

The implication of (3) is that the structural property c(0) (and indeed, 
c(r) for a small range of r where g(r) is itself practically zero) is deter- 
mined almost entirely by therrKJdynamic functions. The generalized mean 
spherical approximatiot. is therefore a model which demands (^^) and (Ic) 
and in consequence sets 
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- -30p/ap)^+ 20U/N) . (8) 

It Is Chen clear chit Che usual mean spherical approximation Is but 
one possible realization of (3). Its range of applicability Is normally 
constrained by Che condition 

3 0p/3P)^ » 2l3U/Nj (9) 

which might be compared with (7), the condition governing Che range of 
applicability of Che generalized mean spherical approximation. Condition 
(9), In turn Is compatible with Chose u(r)'s characterized beyond the hard 
core by either a weak long range tall, or by a stronger tall but with restricted 
range. For the former, the tall has only minor effects on the structure of 
the underlying hard core system and Che application of Che MSA Is then very 
much In Che spirit of thermodynamic perturbation theory^. But the exact 
g(r) reflects, of course, Che entire potential and If the hard core portion 
Is to play any physical role at all, It Is necessary that g(r ■ a"*") — 5, 
where 6 Is significantly different from zero. In fact, of course, the MSA 
leads CO Jump discontinuities In Che solutions for g(r) and c(r) across 

r ■ 7 whose magnitudes depend on a, u(r), and p. We shall need Co Cake this 

Into account In what follows, and shall Introduce an approximation one of 
whose purposes Is to eliminate such discontinuities. 

In proposing ar. extension to the MSA for non hard-sphere systems we 
cannot be guided by the usual arguments founded on (9). Instead any practical 
realizations that might emerge must be sought from Che C^SA postulate as 
summarized by equatloo (3). Such an argument can be given as follows: 
consider a real physical potential possessing not a hard core but a soft 
short range repulsion (examples might be r ”, (ar br and so forth). 


As Is well known, such systems behave at high densities as if their pair 
potentials actually possess an effective hard core of diameter defined 

operationally by g(r • ■ 6 « 1. But though hard-core like . g(r) for 

th«se real systems is certainly devoid of discontinuities, and the same is 
true for c(r). 

The OISA models this real physical problem by caking the u(r) of equation 

(1) to be Just Che portion of the real potential chat lies beyond an assumed 

hard core of dimension J (this 3 often being close, on physical grounds, to 

7 The problem is then solved and as noted, the . pproximate g(r) and 

eff 

c(r) may exhibit jump discontinuities. Let us compare these solutions (which 
obviously depend on the choice of 7 being made) with another approach: 
suppose we return to the tlginal problem and augment the soft potential u(r) 
with a genuine interior hard core whose diameter 7 is chosen so chat ^ ^eff 
Physically (not analytically) such an addition can have Mttle influence either 
in the structural properties (g(r) or c(r), for instance,) or on the high 
density thermodynamic functions. An example of this can be seen in figure 1 
where the results of introducing a hard core into the Coulomb potential, at 
high density, are displayed. Note the small change in c(r), and remember that 
c(0) reflects the thermodynamics. (To appreciate the difference between this 
and when the tail serves as a perturbation see Fig. 29 of ref. 1.) Tnus again 
from the physical viewpoint, we may even choose to regard the results of 
existing Monte Carlo calculations for soft-potentials as in fact equivalent 
to the results of similar calculations, were such calculations actually carried 
out with an Interior hard core (diameter 7 < -g£f) present in the potentials. 
Even though 7 may be quite close to 7^^^ we are evidently approaching a picture 
where the addition of an interior-hard-core is manifested only as a rather 
mild perturbation, a view which is quite contrary to the usual notions governing 
the application of the MSA^. 
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Clearly we can extend this argument and state that any solution of the 
MSA, (equation (1)) which satisfies 

g(r) - 0, r < a , 

g(r ■ a'*’) - 0 , 

(l.e. no discontinuity at r * a) can therefore also represent a MSA result 
for a soft-potential u(r). Such a model can be termed 'a soft-mean spherical 
approxlmut^on' (SMSA) and simply takes the view that for a dense fluid c(r) ■ 

-^u(r) (the simplest realization of (3)) whenever g departs from zero. As 
will become clear later, solutions for g having this character are not guaran- 
teed, but If a solution of SMSA does exist and yields a physically acceptable 
7, this solution should be very similar to the results flowing from the HNC 
approximation (see below). In fact there Is an analytic solution of the SMSA 
for the one component plasma^. For the Yukawa potential an analytic solution 
can also be obtained, the starting point being the solution of the MSA as aug- 

g 

men ted with Yukawa closure . 

In ooth of these cases the SMSA range Is determined by expanding c(r) 
around r ■ 0 In powers of r, and establishing through the choice of 7 where 
Che coefficient of the linear term vanishes. This then guarantees that the 
discontinuity In Che MSA solution for g(r) actually vanishes vl*e* g * 0). 

Since both Coulomb and Yukawa potentials possess Fourier transforms. It Is 
possible to obtain solutions of the SMSA which can be extended to the limit 
of low densities, where 7 -* 0. An Important potential In Che same general 

4 

class Is u(r) erfo(ar)/r whose study by simulation techniques can be of 

some benefit In assessing Monte Carlo results for the OCP via the Ewald image 
9 

method . Analytic solutions of the SMSA for this potential will be particu- 


larly useful. 
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III. MEAN SPHERICAL AND HYPERNETTED CHAIN APPROACHES 

If w« tak« Ilcerally ch« view chat wc have a ayacem wich a hard core 
but Che physical condiclons are such chat Che hard core la not playing a major 
role, Chen we lack the Lmporcanc support of the random phase argvasenc so 
essential to the usual application of the MSA. Ic necessarily follows chat 
any approximations of the mean spherical character must be established by a 
rather different argument. To this end we consider the familiar expansion 
of the total correlation function h(r), given by the diagronsnatic method^^, 
namely 

g(r) - exp (-3u(r) + 9(r) - b(r)] (10) 

where 

9(r) - h(r) - c(r) (11) 

and b(r) (the bridge function) is the negative of the sum of all elementary 
diagrams. The hypemeCCed chain approximation cakes b(r) a 0, so that (10), 
(11) and (2) (the OZ relation) provide together an integral equation that can 
be iteratively solved for a given potential u(r). To make comparison wich the 
MSA, we observe that (10), (11), and (1) also imply a specific choice for the 

4 

bridge function, namely 

( 12 ) 

bj^(r) - h(r) - M g(r) , (r > J) . 

Now in Che SMSA limiting form. 


(r < C) 


g(r) - 0 , 
g(r) ■ 0 , 


(r - a ) 
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which can iclll ba laclaflad by a cholct of b(r) of cha form (12). For lacar 
uaa ho%ravar, we ahall find Ic convanianc to raplaca (12) by cha aingla 
•catemanc 


"SMSA ■ 


(U) 


which formally concinuas Cha definiclon of b into cha raglon whara g(r) 
vanishas. As noted however, this region, by definition of cha model, plays 
no physical role. The utility of (13) lies rather in the following: any 

statement on the behavior of the bridge fxinction is entirely equivalent to 
the specification of an effective interparCicle potential for use in a 
hypemeCCed chain approach. We may refer to this type of presentation of the 

4 

problem as a modified HNC scheme . It is apparent that the SMSA introduced 
above can therefore be cast into a well defined integral equation method. 

For potentials possessing Fourier transforms (the OCP provides a clear 
example) this integral equation has the interesting feature chat it can dis- 
play two classes of solutions that branch from the low density limit. One is 
the "Debye-Hiickel" class characterized by c(r) ■ -9u(r) for all r. The other 
has Che desired high density behavior in g(r), namely a range where g(r) > 0 
which is made self consistently possible by virtue of the correspo iding 
behavior in b(r). A numerical method for finding the desired solution will 
be based on a variational principle for the SMSA (see section V). 

For numerical treatment the SMSA(eq. (13))should still be viewed in the 
context of the limit g(r ■ ) ■ 0 imposed on the modified HNC equation with 

Che bridge function of eq. (12). Observe also that the Percus-Yevick equation 
is characterized within Che modified HNC scheme by the choice 

bpy(r) ■ g(r) - 1 - in g(r) - [c(r) + 3u(r)] , 
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and as is wall known can be cast Into a well defined dlagranmacic expansion. 
The non-Debye-Hiickel branch of the SMSA (eq. (13)) cannot however be cast 
into a diagramnatic expansion. 

Aswlthall Integral equations so far Introduced, It Is not possible to 
say a priori whether any physically ar ceptable solution will In fact emerge 
for an arbitrary potential. But In contrast to other methods, we can show 
that the SMSA solution (provided of course It exists) has very Interesting 
physical features for the class of potentials u(r) whose Fourier transforms 

likewise exist. 

IV. THERMODYNAMIC FUNCTIONS IN THE SOFT .MEAN SPHERICAL APPROXIMATION 

Ue consider a system of particles Interacting via a pair potential that 
Is regular (l.e. lacks a hard-core). If the potential for one additional 
particle is scaled by X (0 < X < 1), and g(r,X) ■ 1 + h(i,X) is the pair 
distribution functlor relative to this particular particle when the potential 
la thus scaled, the excess chemical potential for the system can be *4Tltten^^ 

1 00 

• P ’ dx' dr g(r,X)i3u(r) . (14) 

o o 

Let Y(r,X) be the difference between 3Xu(r) (the potential) and the potential 
of mean force for scaling X. Then 

g(r,X) - exp[-3Xu(r) + Y(r.^)l (15) 

10 

and according to the dlagranraatlc expansion 

y(r;X) » h(r,X) - c(r,X) - b(r,X) (16) 

where, again, b(r,X) Is the negative of the sxim of all elementary graphs. 
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In ch^8« «qu«tlons, w* hav« g(r,X ■ 0) ■ 1, h(r,X • 0) ■ 0, c(r,X ■ 0) ■ 0, 
and b(r,X • 0) » 0. It is also undarstood that g(r,X ■ 1) • g(r), h(r,X ■ 1) 
h(r), and so on. 

We now differentiate (IS) with respect to X and obtain 

3u(r) g(r,X) ■ ■ ^ g(r.X) + g(r,X) ^ Y(r,X) (17) 


which we insert into (14). With the aid of (16) we then find that 


• -p;’c(r)dr p. ^dX’dr h(r,X)r^ (h(r,X) - c(r,X)] 

O 


oX 


(18) 


pJ^dXjdr g(r,X)^ b(i,X) 


which is a simple generalization of 


3u 


ex 


pi'c(r)dr + Pj 


^dXj'dr h(i,X)^(h(r,X) - c(r,X)l 


12 13 

which in turn is a well known ’ expression for the excess chemical poten- 

12 13 

tial in HNC. We can therefore repeat the standard manipulations ' on the 
first three terms of (18) to arrive at 




ex 


- -pj'c(r)dr ^pj'dr h(r) [hlr)-c(r)] - pj* dXjdr g(r)^(r,X)» (19) 


In the HNC approximation (b * 0) the chemical ..otentlal can be calculated 
directly: it is not necessary to integrate the energy equation of state. 
Notice that for any approximation, expression (19) will be equivalent to the 
energy equation of state. 


We now write (19) in a form that is slightly more general and will later 
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permit uf CO make an application to the OCP. We use (3) and (6) (or their 
equivalents for the OCP) and find after a little manipulation that 


au** ■ ^3(|p) ■ ^ au(r)]dr 

+ I" j'h^(r)dr - pj'^dXj'dr g(r,X)^ b(r,X) 

o 


( 20 ) 


which so far remains exact. We next consider the consequences in (20) of 
Che approximations flowing from the statement (13) of the SMSA. (It must 
bo noted that in any approximate theory, aOp/^P)^ (3U/N) in (20) should 

be taken from the compressibility and energy equations of state, respectively.) 

For Che SMSA .e Cake (13) in its scaled form 


bsj^A^*"*^) • h(r,X)- *n g(r,X) (21) 

which is compatible with (3), and insert this into (20). Then using (3) 
we obtain 

'] * f <“> 

a relation which holds however only fer potentials that possess Fourier 

14 -1 

transforms . One such is Che Coulomb interaction (u(r) r ) for which, 

as is well known 



1 

3 



(23) 


Thus tor Che OCP (22) can be rewritten in terms of the Helmholtz free energy 
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a 

N 


SH5»V 




^3P/, 


I + 




(24) 


This lasc txprtsslon demonscracci i unique property of the model: the energy 

equetlon of state and the coapre'islblllty equation of state are completely 
given one In tenaa of thu other. It Is not necessary to know any structural 
details of the solution. Fur':her» If 


a - (3/4 ttP)^^^ 


(25) 


Is the Wlgner-Seltz radius for the plasma, and 
r • (Ze)^3/a 


Is the standard plasma parameter, then In terms of T 




/N). 


(26) 


(27) 


Equations (7), (24) and (27) then constitute a compact and relatively 

simple set of equations connecting the four quantities (3U/N), (3F*^/N), 

3 0p/dP)^ and c(0), and are useful notwithstanding the relative simplicity 

of the analytic expressions for the OCP. 

V. A VARIATIONAL APPROACH TO THE SMSA 

We start with the usual coupling constant expression for the excess free 
15 

energ; 


AW ^ 

3F /N - ^Pidr g(r,X)3u(r) 
• ^ 


(28) 


where g(r,X) is the radial distribution function for Isochorlc systems of 
oartlcles In which the pair potential Is Xu(r). Following the procedure of 
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Morlca and Hirollca^^, and using tquacton (13) (the SMSA statement on the bridge 
function) we obtain 


2P “S’ I‘**itPc(k) + Jn(l-pcOc))] 
~ (2tt) ~ 


(29) 


an approximate form for the free energy equivalent to that obtained frv a 
integrating the energy equation of state. In (29) c(k) is the Fourier 
transform of the direct correlation function c(r). From (28) we can deter- 
mine the excess entropy 


(S®*/Nkg) - -(3F®*/N) + g(r)Su(r)dr . (30) 

But notice that the basic assumption of the generalised mean spherical 
approximation (equation (3)) allows us to write an approximate form for the 
excess entropy 


(S®^Nkg) - -^p^dr h(r)c(r)- ^^^,'dk[pc(k) + .in(l-pc(k))] , 

(2n)'^ 

“ " Jdk[pc(k)H«n(l-oc(k))]j. 


(31) 


We now observe that a functional derivative 

6(S®*/Nkg) 

6c(k) 


of t.ie right hand side of (31) will vanish, provided h and c are connected by 


ph(k) « pc(k)/(l - pc(k)) . 


(32) 
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But (32) is precisely the Fourier cransform of the Omstein-Zernike relation 
(equation (2)). It ther< '*ore follows that though intrinsically approximate 
the entropy functional (31) provides an expression for the excess entropy of 
the system (via the energy equation) that is exact for the SMSA model. Die 
exact solution for the corresponding c(r) is that direct correlation function 
maximizing the entropy functional under the restrictions of the model ; 

Die functional given by the right hand side of (31) is already known^^ 
to possess the property that at the extremum, the Omstein-Zernike relation 
will be satisfied. Diis suggests, in turn, a numerical procedure for obtaining 
a solution to the SMSA which is expected to be valid for any potential. Diis 
procedure is: (i) choose a value for 7 and start with an assumed form for 

c(r), for example 


c(r) ■ c (r) - -3u(r) (r > J) 

09 

c(r) - c^(r) - a + b(l - r/j) -t- (I - r/a) E a P (2r/cT - 1) , 
< n n ’ 

n«0 


(r< 0 ), (35) 

where the are all variational parameters, and the are Legendre poly- 

nomials. Die quantitie-j a and b are fixed by the requirements of continuity 
on c(r) (and its first derivative) at r ■ J . 

(ii) Solve (for a^) the set of equations resulting from the variational 
conditions 


5a\Nk_, 
n B 


0 


(n • 0, I, . . .) 


(ill) alter the value of ? and repeat these steps until 


g(r) s 0 


(r < J) 


to the desired numerical accuracy, and g(r > j) > 0. 
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Observe that since the function c ^( t ) Is expected to be a smooth function 

(and one that is not far from linear) we can expect that in a practical case 

a modest set of 5) will be quite sufficient. The existence of an 

analytic solution of the MSA for the Coulomb potential offers an opportunity 

to check both the validity of this numerical procedure and its utility. 

VI. APPLICATION TO THE ONE COMPONENT PUSMA 

The analytic solution of the SMSA for the OCP is given parametrically^^. 

TT 3 

Let a be the relevant inner hard core diameter, and ri - — pa the associated 
packing fraction. Then we define 

K . - o'- 0 • 

2(i-n) ^ 

Q - - 24ti]^ , (39) 

‘■(1-Tl)^ ■' 

and 

a - K^/24n . (40) 


In terms of these we have the solution for the corresponding equivalei.t 
plasma parameter 


r - 



Further 

^ - -t(i + - •^2)a - (1 + %n)V] . 


and 

c(0) 


(1-n)^ 4(1-ti)2 


(I+^)QK _ 

12r| " bO(’i) 


Finally, 


(41) 


(42) 


(43) 


aOp/3p)^ - 20U/N) - c(0) . 


(44) 
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The energl's given by (42) can be repreaenced exceedingly well by Che 
18 


interpolation 


^ - -0.900Sr + 0.2997T^ + 0.0007 


(45) 


which for r ^ 1 has an accuracy of better chan 0.1%. For large T, on the 
ocher hand, the asymptotic form of the analytic solution gives 


^ - -o.9r + - ... 

N 0 


(46) 


which although only a two term result is actually in remarkable agreement with 
Che results of the complete expansion for T values as low as P ~ 1. It follows 
that if we take Che energy to have the form 


iii 

N 


aP 


+ bP^ + c 


(47) 


thsn the solutions of (7) (24) and (27) are readily obtained: 


3F®'*/N - a: 


P + 2bi + c ifi : + d , 


and 


c(0) ■ - •|bP^ - 2c in P + - 2d - 

+ |bP^ + 2c in P + (■ ^+2d+^^ . 


(48) 


(49) 


The values of a, and b corresponding to (46), Chen give, in particular 


;(0) » -1.2P - + ... 


(50) 


lb 


and 

0Op/3P)^ • -0.6r + ^ . (51) 

In general Che OCP energies given by the SMSA are In gcod '’.^reemenc with 
Che HonCe Carlo resulcs. Buc from (24) we Chen see chac Che requlremencs of 
Chemodynaaic conslsCency musC be badly vlolaced (since we expecc, for large 
r, chac 3 (dp/dp -0.4D. IC is a unique properCy of Che SMSA model Chac 
knowledge of iCs equaCion of scace alone enables us Co draw such a conclusion. 

Ic is noc necessary Co appeal Co any deCails of Che solucion. 

19 

Finally ic Ls insCrucCive Co compare Che SMSA and HNC resulcs for Che 
OCP, parcicularly ac T ^ 1. As a rule, Che cwo soluCions give nearly idencical 
resulcs for nearly all quancities of inceresC, buc parcicularly for (3U/N), 
3(3p/3p)^, and c(o). In comparison wich HNC, however, Che SMSA resulcs are 
slighcly shifced cowards Che MC daca. This observacion suggesCs an inCeresCing 

4 

possibilicy for conscruccing a crucial CesC of a general poinc made in Che 
voncexc of Che modified HNC scheme: Ic has been scaced chac since b(r) encers 

Chie scheme as an effeccive poCenCial, we Chen may expecC chac so long as Cwo 
dlfferenC Cheories diverge one from Che ocher only in a scacemenc governing 
Che long range nacure of b(r). Che consequenC differences in Cheir respeccive 
resulcs musC be qulCe small. Ic follows Chac Che emphasis will focus racher 
nacurally on Che shore range naCure of b(r), especially ics behavior in Che 
region of r corresponding co che flrsC peak of g(r). BuC in jusc chis region do 
we find in Che SMSA a racher weak "poCenCial" (i.e. b(r), since '."'•i hard-core 
is playing no role) Chac indeed shifcs Che HNC (b(r) * 0) resulcs in che 
righc direccion, as we expecC (see Fig. 1). 


VII. CONCLUSION 


A connon approach in Che cheJry of clasaical fluids has been co apply 
perturbation theory from a hard sphere reference system for all potentials that 
possessed a hard core part. It is clear from the discussion in the present 
paper, chat the physical role of the hard core part of the potential has 
not always been interpreted correctly. It is true that the gross structure 
of dense fluids is determined by excluded volume effects and for that matter 
all g(r)'s of dense fluids look very much the same. (The proper context of 
incorporating this first order universality of g(r) is, however, provided by 

4 

Che statement of universality of the bridge functions .) But this need not imply 
chat Che hard core part of the potential will always serve as a good reference 
system. To paraphrase the statement made at the end of ref. 6, a successful 
perturbation theory in the presence of strong Coulomb interactions can be 
obtained only if the hard core part will serve as a perturbation. We see here 
a systematic approach to this situation. 
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Figure Caption 


Fig. I The pair distribution function, g(r)i and the direct correlation 
function c(r), for the OCF at P ■ 70 via various choices of the 
bridge function b(r) as employed In the modified HNC equation. 

1 - dashed line, corresponds to b(r) ■ 0 (l.e. pure HNC). 2 - 
full line corresponds to the SMSA. 3 * dotted line, corresponds 
to; b(r < 1.4) ■ «, b(r > 1.4) ■ 0. 4 - dot-dashed line, corres- 

ponds to b(r) • b„-(r; t \ • 0.4), and actually reproduces quite well 
the Monte Carlo results for the OCP at P ■ 70. 

On the scale of this plot, c(r) for all cases considered above. Is 
indistinguishable from -P/r for r > 1.5. Case 3 features a jump 
discontinuity of about 1 for g(r) and c(r), which Is barely detect- 
able on the scale of c(r). 






